EPJ manuscript No. 

(will be inserted by the editor) 



A new measurement of the structure functions Pll — PttA and 
P LT in virtual Compton scattering at Q 2 = 0.33 (GeV/c) 2 

The Al Collaboration 

P. Janssens 1 , L. Doria 2 , P. Achenbach 2 , C. Ayerbe Gayoso 2 , D. Baumann 2 , J.C. Bernauer 2 , 1.K. Bensafa 3 , R. Bohm 2 , 
D. Bosnar 4 , E. Burtin 5 , N. D'Hose 5 , X. Defay 3 , M. Ding 2 , M.O. Distler 2 , H. Fonvieille 3a , J. Friedrich 2 , J.M. Friedrich 6 , 
G. Laveissiere 3 , M. Makek 4 , J. Marroncle 5 , H. Merkel 2 , U. Miiller 2 , L. Nungesser 2 , B. Pasquini 7 , J. Pochodzalla 2 , 
O. Postavaru 10 , M. Potokar 8 , D. Ryckbosch 1 , S. Sanchez Majos 2 , B.S. Schlimme 2 , M. Seimetz 5 , S. Sirca 8,9 , G Tamas 2 , 
R. Van de Vyver 1 , L. Van Hoorebeke 1 , A. Van Overloop 1 , Th. Walcher 2 , and M. Weinriefer 2 



Department of Subatomic and Radiation Physics, University of Gent, 9000 Gent, Belgium. 

2 Institut fur Kernphysik, Johannes Gutenberg-Universitat, 55099 Mainz, Germany. 

3 LPC, Universite Blaise Pascal, IN2P3, 63177 Aubiere Cedex, France. 

4 Department of Physics, University of Zagreb, 10002 Zagreb, Croatia. 
CEA DAPNIA-SPhN, C.E. Saclay, 91191 Gif-sur-Yvette Cedex, France. 
Physik-Department, Technische Universitat Miinchen, 85748 Garching, Germany. 

Dipartimento di Fisica Nucleare e Feorica, Universita degli Studi di Pavia, and INFN, Sezione di Pavia, Pavia, Italy. 

Jozef Stefan Institute, Ljubljana, Slovenia. 

Dept. of Physics, University of Ljubljana, Slovenia. 

Institute of Space Science, RO 76900, Bucharest-Magurele, Romania. 

Received: . . . / Revised version: . . . 

Abstract Fhe cross section of the ep — > e'p'j reaction has been measured at Q 2 = 0.33 (GeV/c) 2 . Fhe 
experiment was performed using the electron beam of the MAMI accelerator and the standard detector setup 
of the Al Collaboration. The cross section is analyzed using the low-energy theorem for virtual Compton 
scattering, yielding a new determination of the two structure functions Pll — Ptt / £ an d Plt which are linear 
combinations of the generalized polarizabilities of the proton. We find somewhat larger values than in the 
previous investigation at the same Q 2 . Fhis difference, however, is purely due to our more refined analysis 
of the data. Fhe results tend to confirm the non-trivial Q 2 -evolution of the generalized polarizabilities and 
call for more measurements in the low-Q region (< 1 (GeV/c) 2 ). 

PACS. 13.60. Fz Elastic and Compton scattering - 14.20. Dh Protons and neutrons - 25.30.Rw Electropro- 
duction reactions 



1 Introduction 

The internal structure of the proton can be studied using 
virtual Compton scattering (VCS). In this reaction (7* p — > 
7p') a virtual photon 7* scatters off the proton p and a 
real photon 7 is produced. The VCS reaction below the 
pion production threshold allows to measure the gener- 
alized polarizabilities of the proton (GPs) (T]|2]|3]|4). These 
GPs are functions of the photon four-momentum transfer 
squared Q 2 and they describe the polarizability locally 
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inside the proton on a distance scale indicated by Q 2 JH). 
Among the six lowest-order dipole GPs, two are an ex- 
tension of the polarizabilities and /5m obtained in real 
Compton scattering (RCS) 0, quantifying the deforma- 
tion of the charge and magnetization distributions inside 
the proton caused by a static electric or magnetic field, 
respectively. 

VCS is accessed through photon electroproduction (ep 
— > e'p'j) as shown in figured] k, p and q are the three- 
momentum vectors of the incoming electron and proton, 
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and the virtual photon, respectivel}lE k' , p' and q' are 
the three-momentum vectors of the outgoing particles. 
Five independent variables are necessary to define the 
kinematics of the reaction, e.g. the modulus of the virtual 
photon momentum q cm and its polarization parameter £, 
the modulus of the real (outgoing) photon momentum 
q' cm and the polar and azimuthal angles of the real pho- 
ton with respect to the virtual photon direction, # 77cm 
and (p, respectively. The five-fold differential cross sec- 
tion d 5 cr/dk'dQ.] ( idQ.q! cm (with k' = modulus of k') will be 
noted d 5 c. 



CM 




Scattering Plane 



Reaction Plane 



Figure 1. Schematic drawing of the ep — > e'p'7 reaction in the 
center-of-mass of the virtual photon and the target proton. 



2 Experimental determination of the GPs 



The GPs cannot be measured directly. In the physical ob- 
servables (cross sections and asymmetries) they appear 
in specific linear combinations, the structure functions. 
There are six independent GPs, and, by consequence, there 
are six independent structure functions |4j. Only three 
of them, denoted as Pll, Ptt and ^LT/ appear at lead- 
ing order in the low-energy expansion of the unpolar- 
ized cross section. The set of all six GPs can be obtained 
only by double-polarized VCS. Such an experiment was 
performed at MAMI recently for the first time using 
a longitudinally polarized electron beam and measuring 
the recoil proton polarization. The analysis of the dou- 
ble polarization asymmetry |S]|9] will be the subject of 
forthcoming publications. These data can also be used 
for the determination of the unpolarized cross section by 
neglecting the beam and recoil proton polarizations. The 
present paper reports the results of the unpolarized anal- 
ysis from this experiment [8\. The same set of structure 
functions has been determined previously in several ex- 
periments at various values of Q 2 H0lllllrT2l . 

The low-energy theorem (LET) for virtual Compton 
scattering is used to expand the cross section in powers 
of^ m (Elandi)): 



The photon electroproduction reaction contains three 
contributions (see figure^. The reaction is dominated by 
the Bethe-Heitler and Born (BH+B) contributions, where 
the outgoing photon is produced due to bremsstrahlung 
of the electron or proton, respectively. The contribution of 
the BH+B process can be calculated exactly based on the 
proton form factors. The GPs make up the VCS non-Born 
part of the reaction [2\. 





Bethe-Heitler VCS Born VCS non-Born 
Figure 2. Contributions to the ep — > e'p'7 reaction. 



1 All variables defined in the center-of-mass of the virtual 
photon and target proton have an index 'cm'. If no index is 
given the variable is defined in the laboratory system. 



dV = dV 



5^BH+B 



(1) 



where tp is a known phase-space factor, d 5 c BH+B is the 
five-fold differential cross section for the BH+B processes 
and M^ B , which contains the information about the GPs, 
is defined by 



f LL / n 

- — (Pll 



»LT 



»LT 



P TT /£) + Plt • 



(2) 



In this equation e>lt and Cll are known kinematical func- 
tions of (| cm , e, # 77 cm and cp (see e.g. ref. H) ror their com- 
plete definition). At fixed £, two linear combinations of 
structure functions, Pll — ^tt / e and Pu, can be deter- 
mined experimentally. The LET method assumes that, since 
the higher-order terms 0(q' cm 2 ) in eq. fl} are small for 
low q' cm (a condition which holds below the pion produc- 
tion threshold), they can be neglected. The cross section is 
measured in an appropriate kinematical region, i.e. cov- 
ering a range large enough in £>ll and vu, here provided 
by a large coverage in 6jy Cm . Then one forms the quan- 
tity M^ B = (d 5 <r - d 5 cr BH+B ) / (cpq' cm ), which is fitted to a 
linear combination of the two structure functions, as ex- 
pressed by eq. Q. 
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Table 1. Parameters of the spectrometer setup for the VCS90b- 
setting: central momentum, in-plane angle (6) and out-of-plane 
angle (</> o P )- 



Parameter 


Spectrometer A 


Spectrometer B 


Unit 


/^central 


645.4 


539.4 


MeV/c 


e 


38.0 


50.6 


deg 


4>oop 


0.0 


0.0 


deg 



3 Experimental setup and event analysis 

For the present experiment the standard setup of the Al 
Collaboration at MAMI was used [13] together with the 
polarized electron beam and the focal plane proton po- 
larimeter. These two items are not detailed here since they 
play no role in the unpolarized analysis. The beam from 
the MAMI accelerator impinged with an energy of 854.6 
MeV on a liquid hydrogen target. The temperature and 
pressure inside the target cell were constantly monitored 
and the beam charge was measured continuously by a 
Forster probe, allowing to determine the experimental lu- 
minosity £ exp with good precision (well below 1%). The 
mean beam current was about 22 }iA. To prevent local 
boiling of the hydrogen, the beam was deflected with an 
amplitude of a few mm and a frequency of a few kHz. 
The scattered electron and the recoiling proton were de- 
tected in the high resolution spectrometers B and A, re- 
spectively. The setting of the spectrometers is given in ta- 
ble [TJ This setting resulted in the central values q cra = 
600 MeV/c, q' cm = 90 MeV/c, e = 0.645 and <p = 180°, 
and the spectrometer acceptance covered the range [70, 
180]° for 6 llcm . 

Each spectrometer contains two sets of double-planes 
of vertical drift chambers for track reconstruction and two 
scintillator planes for timing and particle identification. 
The gas Cherenkov counter in spectrometer B identifies 
electrons. At the proton side such a Cherenkov detector 
was not present, since the focal-plane polarimeter was 
mounted in spectrometer A. 

The distribution of the coincidence time Tab is shown 
in figure [3] The coincident events are selected in a time 
window of 6 ns and a subtraction of the random coinci- 
dences is performed using the events inside Tab £ [-30, 
-15] ns and Tab £ [15, 30] ns. This correction is small since 
the random coincidences contribute to less than 2% to the 
central peak in fig. [3] (note the logarithmic scale). 

Photon electroproduction events are identified by mis- 
sing-mass reconstruction. The distribution of the square 
of the missing mass in ep — > epX, displayed in fig- 
ure [H shows two peaks corresponding to 7 and 7T° elec- 
troproduction. The 7T° peak is present because the ac- 
ceptance, while centered on q' cm — 90 MeV/c, extends 
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Figure 3. Histogram of the coincidence time Tab for events 
within the analysis cuts. The coincident events inside the cen- 
tral peak (distribution 2) are selected and distribution 1 is used 
for the subtraction of random coincidences. 

20000 



15000 



§10000 

o 

O 



5000 



I I 


I 

□ 1 






□ 2 - 






I I 3 






Ki^i— _ — 



10000 
M£ (MeV 2 /c 4 ) 



20000 



Figure 4. Histogram of the square of the missing mass M^- Dis- 
tribution 1 corresponds to the raw coincidences, while in distri- 
bution 2 only events within the central peak in Tab ar e taken 
into account and the random coincidences are subtracted. In 
distribution 3 also the events in the region of the end caps of 
the target are removed. 

above pion threshold. The separation between the two 
electroproduction processes is excellent (about 30 times 
the peak FWHM). For the calculation of the cross section 
only the events with M| E [-1000,2000] MeV 2 /c 4 are 
used. 

A cut on the target length has been applied to remove 
the events from the interactions of the incoming electrons 
with the end caps of the target, which are much more 
dense than the liquid hydrogen itself. Other cuts are nec- 
essary to select the events inside the desired kinematic 
range: |^ cm - 600MeV/c| < 12MeV/c, |^ m -90MeV/c| 
< 15 MeV/c, |e - 0.645 1 < 0.012 and a range of ±12° 
for the out-of-plane angle of the outgoing photon. After 
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these cuts the signals of the scintillators and Cherenkov 
counters were used to estimate the remaining background 
processes, which were found to contribute to less than 
0.5%. Since this is well below the statistical uncertainty 
of the experiment, no cut was applied on these detector 
signals. The count rates were corrected for the detector ef- 
ficiency. However, this correction was very small and did 
not have any influence on the extracted structure func- 
tions. 

4 Unpolarized cross section and extraction of 
structure functions 

For the determination of the cross section the effective 
solid angle of the detection apparatus is calculated using 
a Monte Carlo simulation [14]. This Monte Carlo takes 
into account the detailed geometry of the apparatus, the 
beam configuration, and all resolution deteriorating ef- 
fects, such as the intrinsic resolution of the detectors, en- 
ergy losses in the materials of the target, etc. The events 
are generated according to the BH+B cross section, which 
is used as an approximation of the real cross section of 
the photon electroproduction reaction (the non-Born con- 
tribution will then be incorporated in an iterative pro- 
cedure). Radiative effects are taken into account as ex- 
plained in 1151 : the acceptance-dependent part is included 
in the simulation and the remaining part, due to virtual 
corrections, is implemented by multiplying the experi- 
mental cross section by a constant factor / cor over the 
complete phase space. This factor equals 0.942 for the 
kinematics of the present experiment. The simulation re- 
produces the radiative tail very well, as can be seen in fig- 
ure [5] In this figure the simulated distribution is normal- 
ized using the factor £ exp / (£ s im/cor)/ where £ s im is the 
luminosity corresponding to the simulated events. The 
stability of the experimental cross section versus differ- 
ent cuts is better than 1%. 

The central momenta of the spectrometers were cali- 
brated in absolute using the -distribution. By simul- 
taneously adjusting the experimental peak position (on 
the simulated one) and minimizing its width, one obtains 
the two central momenta. This adjustment was done on a 
kinematical phase space as wide as possible. For the sub- 
set of events within the analysis cuts, it results in a slight 
offset between the experimental and simulated peak po- 
sitions (fig. [5]l, reflecting the uncertainty of the calibra- 
tion. This uncertainty is estimated to be ±3 • 10~ 4 (in rel- 
ative) of the central momentum of each spectrometer. It 
can fully account for the observed offset, i.e. peak posi- 
tions in fig. [5] would coincide by changing either one mo- 
mentum or the other within the quoted uncertainty. 



T 




M£ (MeV 2 /c 4 ) 

Figure 5. Experimental (histogram) and simulated (full line) 
histogram of M^. Both distributions are obtained after apply- 
ing all cuts. The slight offset in position between both maxima 
is discussed in the text. 

We now explain the iteration method used to obtain 
the experimental (ep — > e'p'j) cross section. It is impor- 
tant to use the most realistic cross section for the event 
generation in the simulation fl4l . In a first step, using 
dc BH+B in the simulation, one determines the experimen- 
tal cross section and extracts the structure functions Pll — 
Ptt / £ and Plt, as explained in section [2] Then in a sec- 
ond step the cross section in the simulation is modified 
to include the GP effect measured in the first step. The 
effective solid angle is recalculated, and the whole pro- 
cedure is iterated several times. After three iterations a 
stable result is obtained. The effect of the iterations on 
the solid angle is shown in figure [6]a. The effect is small 
(< 2%), but, due to its pronounced # 77cm -dependence, 
it has a substantial influence on the obtained structure 
functions (see section |5|. Another feature is that the it- 
eration procedure yields the same result, independent of 
the initial value for the structure functions. This is rep- 
resented in figure [6]b. When the BH+B cross section is 
used in the first step of the simulation, the starting point 
is (P LL - Ptt/£, Plt) = (0,0) GeV" 2 . 

The final cross section is displayed in figure and the 
cross section values are given in table|2] The main sources 
of systematic uncertainties are the calibration of the mo- 
menta and angles of the reconstructed particles, the nor- 
malization of the cross section (luminosity), the radia- 
tive corrections and the simulation of the solid angle. To 
study the first point, the central momenta of the spec- 
trometers were changed by ±3 • 10~ 4 in relative and the 
data were re-analyzed, yielding new values for the cross 
section (and the structure functions). The maximal devi- 
ation w.r.t. the original values was taken as the uncer- 
tainty due to momentum calibration. A procedure along 
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Figure 6. a: Influence of the iteration procedure on the effective 
solid angle. The points at left of (0yj Cm = 180°) correspond to 
f = 0°, the points at right to <p = 180°. b: Convergence pat- 
tern for three different starting points, represented by the val- 
ues of ( Pll — Ptt/ £ / ^lt) m parenthesis. The cross at the final 
point gives the size of the statistical uncertainty. The figure is 
obtained with the proton form factors of ref. |16| (see section[5). 



the same lines allows to estimate the systematic error due 
to the uncertainty in horizontal angle (spectrometer angle 
plus transfer matrix), taken equal to ±0.1 mr in each arm. 
These sources of error are $<yy Cm -dependent, changing the 
shape of the cross section. The three other sources are 
#77cm-independent; summed quadraticaHy, they cause an 
error in the overall normalization of the cross section of ± 
2%. The statistical error on the cross section is generally 
smaller, about ± 1.4% for most data points (tableEJ. 

The two structure functions Pll — Ptt / £ and Plt are 
extracted by a linear fit of the quantity Mq 16 / £>lt as a 
function of I'll/ ^lt- They are determined at a fixed value 
of cjcm = 600 MeV/c, or equivalently at a fixed value of 
Q 2 = 0.33 (GeV/c) 2 . The fit is performed via a x 2 mini- 
mization, which also provides the statistical error on the 
structure functions. The result is shown in figure [8] The 
reduced x 2 of the fit, 2.6 for 20 d.o.f., suggests that the 



Table 2. The ep — > e'p'j cross section, in pb/(MeV.sr 2 ), 
obtained in this experiment (after iterations) at fixed q' cm = 
90 MeV/c, q cm = 600 MeV/c and £ = 0.645. The error is statis- 
tical only. The BH+B cross section is calculated using the form 
factors of reference fl6l . 



^77cm 


<P 


d c + 


d 5 cr 


177.5° 


180° 


0.129 


0.146 ± 0.002 


172.5° 


180° 


0.132 


0.137 ± 0.002 


167.5° 


180° 


0.136 


0.140 ± 0.002 


162.5° 


180° 


0.142 


0.148 ± 0.002 


157.5° 


180° 


0.148 


0.150 ± 0.002 


152.5° 


180° 


0.153 


0.155 ± 0.002 


147.5° 


180° 


0.156 


0.154 ± 0.002 


142.5° 


180° 


0.158 


0.156 ± 0.002 


137.5° 


180° 


0.158 


0.159 ± 0.002 


132.5° 


180° 


0.157 


0.157 ± 0.002 


127.5° 


180° 


0.155 


0.154 ± 0.002 


122.5° 


180° 


0.151 


0.142 ± 0.002 


117.5° 


180° 


0.147 


0.140 ± 0.002 


112.5° 


180° 


0.142 


0.135 ± 0.002 


107.5° 


180° 


0.137 


0.128 ± 0.002 


102.5° 


180° 


0.132 


0.120 ± 0.003 


97.5° 


180° 


0.126 


0.119 ± 0.003 


92.5° 


180° 


0.122 


0.109 ± 0.003 


87.5° 


180° 


0.117 


0.102 ± 0.003 


82.5° 


180° 


0.113 


0.103 ± 0.003 


77.5° 


180° 


0.109 


0.098 ± 0.003 


72.5° 


180° 


0.106 


0.099 ± 0.005 


177.5° 


0° 


0.132 


0.149 ± 0.003 


172.5° 


0'- 


0.141 


0.167 ± 0.004 



Table 3. Estimation of the systematic error on the structure 
functions (in GeV~ 2 ) . 





^LL" 


P TT /£ 


Plt 


Momentum calibration 


± 


2.7 


± 1.0 


Horizontal angles 


± 


1.2 


± 0.4 


Normalization of the cross section 


± 


0.6 


± 1.9 


Total systematic error (quadr.sum) 


± 


3.0 


±2.2 



higher-order terms 0(q' cm 2 ) in eq. ((TJ, although small, are 
not completely negligible. This effect of the LET trunca- 
tion has not been considered in the error budget. The sys- 
tematic errors on the structure functions are summarized 
in table |3j where, as for the cross section, they are sepa- 
rated into an angle-dependent and an angle-independent 
part (first two lines and third line, respectively). 
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Table 5. Results for Pll — Ptt/ £ and Pll after applying the iter- 
ation procedure for q cm = 600 MeV/c using different form fac- 
tor parameterizations. The first error is statistical, the second 
the systematic one. 





Pll 


Form 


(GeV- 2 ) 


(GeV- 2 ) 


factor 


This work 27.1 ± 1.9 ± 3.0 


-8.0 ± 0.7 ± 2.2 


m 


28.5 ± 1.9 ± 3.0 


-5.2 ± 0.7 ± 2.2 




28.6 ± 1.9 ± 3.0 


-10.1 ± 0.7 ± 2.2 


m 


27.4 ± 1.9 ± 3.0 


-6.8 ± 0.7 ± 2.2 


1201&12H 



Figure 7. The unpolarized ep — > e'p'j cross section measured 
at Icra = 90 MeV/c, q cm = 600 MeV/c, e = 0.645 and cp = 180°. 
The dotted line shows the BH+B cross section, calculated with 
the proton form factors of ref. fl6l - The full line includes the 
effect of the structure functions obtained in the experiment (first 
line of table [5). The two most left points correspond to cp = 0°. 

Table 4. Results for P L l — Plt/ £ an d Plt without the iterations 
for q cm = 600 MeV/ c using different form factor parameteriza- 
tions. The first error is statistical, the second the systematic one. 
£ in ref. ITTI (e = 0.620) was slightly different from the present 
experiment (£ = 0.645). 



> 

CD 

O 






Pll ~ Ptt/£ 


Plt 


Form 




(GeV" 2 ) 


(GeV" 2 ) 


factor 


This work 


23.3 ± 1.9 ± 3.0 


-6.6 ± 0.7 ± 2.2 


m 




24.3 ± 1.9 ± 3.0 


-3.9 ± 0.7 ± 2.2 


tl7l&IT8l 




24.7 ± 1.9 ± 3.0 


-8.9 ± 0.7 ± 2.2 


ed 




23.7 ± 1.9 ± 3.0 


-5.3 ± 0.7 ± 2.2 


1201&12T1 


Ref. HQ 


23.7 ± 2.2 ± 4.3 


-5.0 ± 0.8 ± 1.8 


|20| 



5 Discussion 

The effect of the GPs in the photon electroproduction cross 
section is small; in the kinematics of the experiment it 
reaches at maximum ~ ±10% (see fig. [7}. Therefore any 
small change at the cross section level may induce a rela- 
tively large change in the fitted structure functions. 

A first example is provided by the iterative calcula- 
tion of the solid angle, explained in section [4] Using this 
procedure, Pll — Pit/ 6 is increased by 16% and Plt by 
20-30% (see tablesHandE). In the first VCS experiment at 
MAMI ITfl , this iterative procedure was not pushed to its 
convergence point, because its effect was smaller than the 
statistical uncertainty, at the cross section level. Indeed, 
as can be seen in figurelH-a, the iterations change the solid 



Figure 8. The linear fit of the structure functions: Pll — Ptt / e 
(slope) and Pll (intercept). Error bars are statistical only. The 
angular values reported along the schematic dotted line show 
the evolution of the points with 8 77cm . The two extreme points 
(677cm, f) = (72.5°, 180°) and (172.5°, 0°) are removed from 
the fit. The figure is obtained with the proton form factors of 
ref. (161 and after final iteration. 

angle (and hence the cross section) by less than 1% in the 
main part of the phase space, whereas the statistical un- 
certainty on the cross section was ~ 2-3%. Therefore the 
result of [11 J is non-iterated. It can be compared to the 
non-iterated result of the present analysis, performed at 
the same value of q cm . As shown in table SJ at this level 
the agreement is strikingly good between the two exper- 
iments. 

A second example is provided by the form factor pa- 
rameterizations (for recent reviews on nucleon form fac- 
tors, see refs. [22] and [23 1). The obtained structure func- 
tions depend on the choice made for the proton form fac- 
tors Gg and Gj^, because these observables enter the cal- 
culation of the BH+B cross section. In this analysis vari- 
ous form factor parameterizations have been considered 
|[T6Tl7inillT9ll2TJll2Tl (see figure O. Mainly P LT is sensitive 
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to this choice: going from the parameterization of ref . [20 J 
(or (22) to the other ones, Plt changes by (^43)% while 
Pll — Ptt/ £ changes by 1-4% only (cf. table 0. This is 
caused by the fact that a change of form factor induces 
mainly a variation of the magnitude of the BH+B cross 
section without affecting too much its # 77cm -dependence 
(see figure[9}, resulting mostly in a change of the intercept 
in figure [8] 
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Figure 9. The BH+B cross section, at the kinematics of the 
experiment, for proton form factors as parameterized in |16| 
(dotted), HZJ-CHI (solid), QU (short-dashed) and (20) (long- 
dashed). The corresponding form factor parameterizations 
are shown in insert, with G D = [1 + (Q 2 /0.71 (GeV/c) 2 )]~ 2 . 
Ref. l2Tl gives curves very similar to 1201 . 

It should be noted that in flT] the experimental cross 
section was compared to the theoretical one also at very 
low q' cm . This test favored the form factor parameteriza- 
tion of ref. E0l which was consequently chosen in the 
analysis, and an overall form factor uncertainty was em- 
bedded in the systematic error^. The present experiment 
is only performed at q' crn = 90 MeV/c, preventing such 
normalization test at low q' cm . The uncertainty due to the 
proton form factors is not included in the numerical value 
of the systematic error of table [5] It is represented ex- 
licitely by the four different lines of results in this table 



The resulting structure functions are displayed in fig- 
ure [10] together with the other existing measurements. 



2 This uncertainty is accounted for in the last line of table [4] 

3 If however one wants a single number for this system- 
atic error, one may take the half-difference between the two 
extreme results of table [5) i.e. ±0.7 (resp.±2.4) GeV~ 2 for 
the first (resp. second) structure function. This would give af- 
ter quadratic sum a total systematic error of ±3.1 (resp. ±3.3) 
GeV~ 2 on P LL - Pjj/e (resp. P LT ). 
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Figure 10. The structure functions measured in this experi- 
ment (•) (with the form factors of ref. |20|) and previously at 
Bates (TO), MAMI QJJ, JLab |TJ) and in RCS (6). The solid curve 
is the HBChPT calculation (24 j. The dashed curve is an example 
of the dispersive (DR) calculation 1251 with dipole ansatz pa- 
rameters (A a , Ap) = (1.8, 0.5) GeV. Curves are calculated for e 
= 0.645. Data points correspond to different £ (0.90 (TO), 0.62 [TTJ, 
0.645 (this exp.), 0.95 and 0.88 (TJ)). They are obtained from ei- 
ther a LET analysis or a dispersive one (DR). Some points are 
shifted in abscissa for visibility. The inner error bar is statistical, 
the outer one is the quadratic sum of statistical and systematic 
errors. 

At low Q 2 they can be compared to a calculation in the 
framework of the heavy baryon chiral perturbation the- 
ory (HBChPT) ||24|| . The calculation at order p 3 for the 
structure functions evaluated at Q 2 = 0.33 (GeV/ c) 2 and 
£ = 0.645 gives P LL - Prr/e = 25.5 GeV" 2 and P LT = 
—5.3 GeV" 2 , in overall good agreement with the values 
measured in this experiment. At next order in HBChPT, 
only the spin GPs have been evaluated I26II27I . A com- 
plete ChPT calculation of the present structure functions, 
combining scalar and spin GPs, is still to come. The dis- 
persive formalism of refs. [25.28] offers an alternative ap- 
proach in which the electric and magnetic GPs can be fit- 
ted from the experiment. For a more extensive compari- 
son to the theoretical models we refer the reader to e.g. 
ref. El. 

In the structure functions presented here, the contri- 
bution of the spin GPs (and hence Pjj, which contains 
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only spin GPs) is expected to be small. Therefore the data 
for Pll — Ptt/ £ essentially reflect the behavior of Pll, 
which is proportional to the product Gg ■ ag. The data 
do not follow a simple dipole shape over the full mea- 
sured Q 2 -range. To fix the shape, more measurements are 
needed in the region up to 1 (GeV/c) 2 . For the structure 
function Plt, the scalar part is proportional to Gg • /3m- 
The measured points tend to confirm the existence of an 
extremum of /3m at low Q 2 , traditionally explained by the 
interplay between diamagnetic and paramagnetic contri- 
butions. To summarize, both structure functions show a 
non-trivial Q 2 -behavior, that can be related to the pion 
cloud structure of the nucleon. However, experimental 
data are still scarce and more measurements in the Q 2 - 
region between and 1 (GeV/ c) 2 would help to gain in- 
sight into the matter. 

6 Conclusion 

New values for the structure functions Pll — ^tt / e an d 
Plt were obtained from the present VCS experiment per- 
formed at MAMI. Apart from a well-understood system- 
atic difference, the results are in very good agreement 
with the ones obtained in the previous MAMI experi- 
ment ITlTI . They are also in good agreement with the cal- 
culation of HBChPT |24| . The new feature in the present 
analysis is an iterative procedure in the calculation of the 
solid angle, inducing an increase of both structure func- 
tions. The effect of different proton form factor param- 
eterizations has also been studied in detail. More pre- 
cise measurements of these form factors at low Q 2 ||30l 
[311 will help to reduce the uncertainties in the determi- 
nation of the GPs. More VCS measurements at low Q 2 
would help to investigate the non-trivial behavior of the 
GPs and their connection to the mesonic structure of the 
nucleon. 
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